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. Abstract. We analytically calculate the average value of i-th largest Schmidt 

^ ' coefficient for random pure quantum states. Schmidt coefficients, i.e., eigenvalues of 



the reduced density matrix, are expressed in the limit of large Hilbert space size and 
O^l , for arbitrary bipartite splitting as an implicit function of index i. 
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\0 • Entanglement is a valuable quantum resource and as such a lot of attention 

^ ■ has been devoted to studies of entanglement and ways of exploiting entanglement 



present in various quantum states. Random quantum states are interesting for several 
reasons. They are used in various quantum information processes like e.g. in quantum 
communication. Also, a generic quantum evolution falling under the label of quantum 



chaos will produce a random state from almost arbitrary initial state. It is therefore 
important to understand entanglement properties of random quantum states. It comes 



as no surprise that the subject has a long history, for instance, entropy of a subsystem has 
been calculated in [H El EJ HI [5] while purity has been studied in [H [6] . From a practical 
point of view a relevant question is how to generate random pure states? Several methods 
have been proposed, using the so-called pseudo-random unitary operators [3, El El [ID] 
which are characterized by the number of random parameters scaling polynomially with 
the number of qubits and not exponentially as it would be required in order to generate 
true random unitary operators. Of course, using less than a full set of exponentially 
many parameters causes pseudo-random operators to be just an approximation to the 
true ones. Nevertheless, if a common criterion of entanglement such as purity or von 
Neumann entropy of a subsystem is used, the convergence to the asymptotic values of 
purity or entropy is observed already after polynomial number of steps. While purity 
might be a perfectly good benchmark for certain purposes we must nevertheless keep in 
mind that purity or entropy of a subsystem is just one of many possible entanglement 
measures of pure states. In fact, to fully characterize entanglement properties of a given 
bipartite state a full list of Schmidt coefficients is needed. If we have a bipartite Hilbert 
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space H = Ha ® Hb of dimension NK, with being dimension of Ha and K of Hb, 
the Schmidt decomposition is 

7V-1 , 



where and \wf) are orthonormal on respective subspaces Ha and Hb- Squares 
of Schmidt coefficients Aj are nothing but the eigenvalues of the reduced density 
matrix, p = tr ?/')(?/' |, obtained by tracing over subspace Hb- Without sacrificing 
generahty we will assume N < K as well as that the eigenvalues Aj,z = 0, . . . , N — 1 
of p are listed in decreasing order. Schmidt coefficients are invariant under local 
transformation (transformations acting separably on Ha and Hb) and fully specify 
bipartite entanglement of state {ip)- Therefore, a set of — 1 entanglement monotones, 
for instance, sums of k smallest eigenvalues, is needed to characterize pure state 
entanglement [11]. A sum of k smallest eigenvalues is also a good estimate for the 
truncation error in recently proposed numerical method [12] for simulation of weakly 
entangled quantum systems using matrix product state representation [13]. For random 
states joint distribution Pn,k(Ao, • • • , ^n-i) = -Pn.kI-^) of eigenvalues is a well known 
quantity [21 [6], see equation (JSj) below. However, it is not particularly convenient 
to deal directly with the distribution in A^ dimensional space. It would be much 
simpler to have some scalar quantities which could be used as a benchmark for the 
randomness of quantum states. Because the eigenvalues fully determine entanglement 
we are going to calculate the average value of Aj for random quantum states, i.e., the 
average value of i-th largest eigenvalue of the reduced density matrix, (Aj), where the 
angle brackets denote averaging over random states. Random vectors in H can be 
obtained in two equivalent ways. As columns of random unitary matrices distributed 
according to the invariant Haar measure, which can in turn be obtained by, e.g., Hurwitz 
parametrization of the unitary group. A second, simpler procedure, is to simply draw 
real and imaginary parts of expansion coefficients of as random Gaussian variables, 
subsequently normalizing the vector [6]. We will be especially interested in the limit of 
large dimensions A^. Because Aj scale as ~ 1/A^ (due to = 1); we will define a 

function /(x), such that 



In the limit N ^ (yo the function f[x) becomes a continuous function with the definition 
range x G [0, ij^l- Deriving and studying the function f{x) is the main aim of this paper. 
There are several ways of calculating (Aj). The most direct one is by using joint density 
of eigenvalues PAr,i<-(A). While such method is exact for all A^ it is rather complicated. 
As we are mainly interested in the asymptotic A^ — > oo behaviour we are going to use 
simpler method to derive f{x) in the second part of the paper. 

I The reason to set x = {i + l/2)/N and not e.g. x = i/N is purely empirical. By comparing our 
analytical solution for f{x) with numerics for finite N we have found that setting x — {i + l/2)/N gives 
better agreement than e.g. x = i/N or x = (i + 1)/N. Note though, that the difference between all 
these choices of x goes to zero with N oo. 




(1) 




X = 



i + 1/2 
N 



(2) 
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For random \ip) the distribution of eigenvalues is [21 E] 

To calculate the average of i—th largest eigenvalue Aj or its distribution it is useful to 
first define the probability E{i;s) that exactly i eigenvalues lie in the interval [0, s]. 
Then the probability density of i-th eigenvalue is given by s), obtained as 

F(0; .) = F{^■, s) - F(z - 1; .) = for ^ > 0. (4) 

For standard Gaussian random matrix ensembles the distribution of the largest 
eigenvalue has been calculated in [HI [15]. Average value of i-th eigenvalue is then 
simply (Aj) = JsF{i; s)ds. For instance, the smallest and the largest eigenvalues are 
given by 

(A„,in) =N f'dss [' PN,Kis,\)dX 

Jo Js 

(A^ax) = N f'ds s r PnM^, A)dA, (5) 
Jo Jo 

where the inner integral is over N — 1 arguments of Pn,k ©• Using these formulae 
in the simplest case of a symmetric bipartite cut with N = K one obtains for N = 2 
(Ao) = I and (Ai) = |, while for = 3 one gets (Aq) = fg, (Ai) = and (Aa) = ^. 
We have calculated (Amin) for several other small values of N with the result always 
being We therefore conjecture that for N = K we have in general 

(Amin) = (6) 

Even though the integrals in ([5]) can be connected to Laguerre polynomials [16], [6] they 
are non-trivial because the inner integrations are over a finite range. To calculate f{x) 
in the limit N ^ oo we will use a different approach. Instead of using the exact Pn,k 
and afterwards taking the limit we will instead use the density of eigenvalues p{X) which 
is already evaluated in the N, K —>■ oo limit. 

Eigenvalue density p{r) of scaled eigenvalues r = AjA^ can be obtained by 
integrating out A^ — 1 eigenvalues in Pn,k- For large A^ and K the result is [3] 

Pi-) = ' (7) 

if r G [a, b], otherwise p{r) = 0. Parameters a, b are given as 

a=(l-v^)^ 6=(l + v/^)^ w=^. (8) 

The above density of eigenvalues ([7]) is well known in statistics, where it occurs 
as eigenvalue density of random covariance matrices WW"^ , also known as Wishart 
matrices, and is called Marcenko-Pastur law [17]. In the simplest case of an equal 
partition, N = K, we have a = and 6 = 4, resulting in p(r) = ^\J'^/t — 1- Moments 
of eigenvalue distribution, (A'^) = tr p"^ have been calculated in ^18j. The calculation 
of f{x) from p(r) is simple because f{x) is by definition a decreasing function of x. 
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Noting that p{t) = l/\f'{x*)\, where x* is such that f{x*) = r, we can immediately 
write differential equation for f{x), 

fix) = -1/pix), (9) 

with f'{x) = df{x)/dx and the boundary conditions /(O) = b and /(I) = a. Let us first 
study the symmetric case of equal bipartition, N = K, and then separately the general 
N < K. 

Integrating Qj once, using ([7]) and setting N = K, we get f{x) in an implicit form 
/(a;)=4cosV, |x = - ^ sin (2(^), {N = K}. (10) 




Figure 1. Theoretical f{x) = NX{x) (fTO| for N = K. In tlie riglit frame is a semi-log 
plot. 



Functional form of f{x) (JTUj) can be seen in figure [TJ For small or for large values 
of X, i.e., for large and small eigenvalues, one can expand implicit solution ffTOl) to 
get f{x) ^ 4 - 4(37rx/4)^/'' for x < 1 and f{x) ^ 7r^(l - x)^/4 for 1 - x < 1. 
We see that the largest eigenvalue is (Aq) ~ 4/A^, while the smallest one scales as 
(Aat-i) ~ in accordance with the conjecture ([6]). Theoretical formula for /(x) 

ffTOj) gives asymptotic N oo result. What about the eigenvalues for finite A^? We 
performed numerical calculation of average Aj by randomly drawing vectors and then 
calculating the eigenvalues of the reduced density matrix. The convergence to the 
asymptotic f{x) is pretty fast and already for moderate iV = 2" with n = 4 {NK = 2^), 
theoretical formula (ITOll agrees with numerical Aj to better than 1%. In figure [2] we 
can see that the relative error decreases as ~ apart from the smallest eigenvalue 
(and possibly few smallest ones) for which theoretical formula ( fTOl) gives Xn~i ~ 0.62/A^^ 
while the exact value is conjectured to be ([6]). We have also numerically calculated 
the width 6\i of the distribution of Aj, ((5Aj) = (A?) — (Aj)^. Numerical data (not 
shown) indicate that the width 6Xi is approximately equal to one half of the spacing 
Aj — Aj_i. Using theoretical formula for the spacing (i.e. for f'{x)) we have found that 
an approximate relation 6Xi ~ A/ [N"^ ^ A/ f[x) — 1) holds, where /(x) is given in (ITOl) . 
With increasing the relative width decreases and therefore the eigenvalues become 
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Figure 2. Relative error of numerically calculated average Xi for N — 2", n = 4, 5, 7, 9 
(top to bottom), as compared to the theoretical formula for f{x) pOjl . Relative error 
decreases as 




Figure 3. Theoretical error ri{x) (I12I10|) committed by retaining a fraction x of the 
largest eigenvalues. In the inset is a semi-log plot of the same ri{x). All for N = K. 



increasingly peaked around their average (A,). If we want to approximate a given state 
with less than a full set of parameters, e.g. as a matrix product state [I3], the quantity 
that estimates an error made by retaining only a fraction x of largest eigenvalues is ri{x), 
defined as 

ri{x) = dx X{x). (11) 

Jx 

Changing integrating variable from x to (p the above integral can be evaluated. Using 
our f{x) we get 

r/(x) = l--{2<^-lsin(4<^)}, (12) 

with ip being the solution of equation (fTOj) . Plot of ri{x) (fT2l) is shown in figure [31 An 
interesting thing we can see is that because f{x) behaves as ~ (x — 1)^ close to x = 1, 
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the smallest eigenvalues contribute very little to r]{x). For instance, the first 80% of 
eigenvalues sum to more than 0.993, i.e., discarding 20% of the smallest eigenvalues will 
cause fidelity error of no more than 1 — F ^ 0.007. 

Let us now go to the general case of N < K. By a similar procedure as for 
symmetric N = K situation we arrive at the implicit formula for f{x), 

f{x) = a + {b-a) cos^ ip, {N < K} 

71 1 + w 1 1 — w {fa \ 

-X = - ^ sm (2^) - ^ arctan (^^ - tan j , (13) 

where a, b and w are given in ([H]). By comparing matrix eigenvalue equation for Wishart 
matrices (i.e., W^W^^, which is nothing but the unnormalized reduced density matrix in 
our language) and zeros of Laguerre polynomials in the limit N,K ^ oo, it has been 
shown in [21] that the two converge to one another. Therefore, in the limit N, K ^ oo 
the i-th zero of a generalized Laguerre polynomial ^'^^\Kx) converges to NXi, i.e., 
to f{xi = {i + l/2)/N) with /(x) given in equation (fT3l) . 




Figure 4. Scaled eigenvalues \{x) for a bipartite cut to n/2 — r and n/2 + r qubits. 
Circles are numerics for n ~ 14 while the full line is theoretical f{x) (fTS]) . 

To verify theoretical formula f|T3l) we performed numerical calculation for random 
vectors in 2" dimensional Hilbert space, making a bipartite cut to n/2 — r and n/2 + r 
qubits, i.e., = 2"/^"'', K = 2"/^+*" and w = 1/2^''. Comparison of numerics and theory 
for r = 0, 1, 2, 3 and 4 is in figure HI One can see that with increasing r, as the subspaces 
TYa and Tie become of increasingly different size, the reduced density matrix on TYa has 
increasingly flat spectrum. When one of the subspaces is much larger than the other, 
that is when w = N/K <^ 1, we can simplify the general equation for f{x) ( |T3i) to arrive 
at 

/iV TT 1 

/(x)^l + 2y-cos(2<^), -x = (p--sin(4(^) {N <^ K}. (14) 

For small w the eigenvalues of p therefore deviate from 1/A^ by ~ jfy/w. The fact that 
with increasing r all eigenvalues are essentially ~ 1/A^ is well known. Namely, if we 
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trace a random pure state over large "environment" the resulting density matrix will be 
almost proportional to the identity, e.g., its purity is close to [H |19]. In addition, 
deviations from this totally mixed state go to zero if K increases at a fixed A^. This is a 
consequence of the so-called measure concentration in many dimensional spaces [20j 




Figure 5. Theoretical error r]{x) (|12ll3p for the same Hilbert space sphtting as in 
figure m In the right frame rj is plotted as a function of y = xN/ \/NK so that a given 
y represents the same number of eigenvalues for all r and not the same fraction, as in 
the left plot. 



The error r]{x) (fTTj) made by retaining only xN largest eigenvalues is for N < K 
given by exactly the same expression as for N = K fll2l) . the only difference is that ip is 
given by the solution of f|T3|) . Because the spectrum of p becomes flat for large r, rj{x) 
approaches linear function. This can be seen in figure [5l As the number of retained 
eigenvalues xN = x2"/^~'" depends on r, we also plot in the right frame of figure [5] 
dependence of rj on the absolute number of retained eigenvalues. From this plot one can 
see that for instance, keeping 0.2y/ NK largest eigenvalues in all bipartite cuts, we make 
an error rj ^ 0.45 for symmetric cut with r = 0, error rj ^ 0.40 for r = 1, r/ ^ 0.12 for 
r = 2, while for larger r's the error is zero because the dimension of TYa is smaller 
than the number of retained eigenvalues {i.2\/NK. 
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§ For uniform measues in high dimensional space deviations from the average values are very small. 
As a simple example: probability that the sum of d bounded random numbers deviates by more than 
e from its average is exponentially small in the dimension d as well as in e^. 



Entanglement of random vectors 



8 



References 

[I] Elihu Lubnik, Entropy of an n-system from its correlation with a k-reservoir, J. Math. Phys. 19, 

1028-1031 (1978). 

[2] Seth Lloyd and Heinz Pagels, Complexity as thermodynamic depth, Ann. Phys. (NY) 188, 186-213 
(1988). 

[3] Don Page, Average entropy of a subsystem, Phys. Rev. Lett. 71, 1291 (1993). 

[4] Jorge Sanchez-Ruiz, Simple proof of Page 's conjecture on the average entropy of a subsystem, 
Phys. Rev. E 52, 5653-5655 (1995). 

[5] Siddhartha Sen, Average Entropy of a Quantum Subsystem, Phys. Rev. Lett. 77, 1 (1996). 

[6] Karol Zyczkowski and Hans-Jiirgen Sommers, Induced measures in the space of mixed quantum 
states, J. Phys. A 34, 7111-7125 (2001). 

[7] Joseph Emerson, Yaakov S. Weinstein, Marcos Saraceno, Seth Lloyd and David G. Cory, Pseudo- 
Random Unitary Operators for Quantum Information Processing, Science 302, 2098-2100 (2003). 

[8] Yaakov S. Weinstein and C. Stephen Hellberg, Entanglement Generation of Nearly Random 
Operators, Phys. Rev. Lett. 95, 030501 (2005). 

[9] Joseph Emerson, Etera Livinc and Seth Lloyd, Convergence conditions for random quantum circuits, 
Phys. Rev. A 72, 060302(R) (2005). 

[10] Roberto Oliveira, Oscar C. O. Dahlsten and Martin B. Plenio, Efficient Generation of Generic 
Entanglement, preprint 'quant-ph/0605126 (2006). 

[II] Guifre Vidal, Entanglement monotones, J. Mod. Opt. 47, 355-376 (2000). 

[12] Guifre Vidal, Efficient Classical Simulation of Slightly Entangled Quantum Computations, 
Phys. Rev. Lett. 91, 147902 (2003). 

[13] Tomaz Prosen and Marko Znidaric, Is efficiency of classical simulations of quantum dynamics 
related to integrability?, preprint |quant-ph/0608057l (2006). 

[14] Craig A. Tracy and Harold Widom, Level-spacing distributions and the Airy kernel, 
Commun. Math. Phys. 159, 151-174 (1994). 

[15] Craig A. Tracy and Harold Widom, On orthogonal and symplectic ensembles, Com- 
mun. Math. Phys. 177, 727-754 (1996). 

[16] Madan Lai Mehta, Random Matrices, 2nd. ed. (Academic Press, 1991). 

[17] Vladimir A. Marcenko and Leonid A. Pastur, Distribution of eigenvalues of some sets of random 

matrices. Math. USSR-Sb 1, 457-483 (1967). 
[18] Hans-Jiirgen Sommers and Karol Zyczkowski, Statistical properties of random density matrices, 

J. Phys. A 37, 8457-8466 (2004). 
[19] Vivien M. Kendon, Karol Zyczkowski and William J. Munro, Bounds on entanglement in qudit 

subsystems, Phys. Rev. A 66, 062310 (2002). 
[20] Patrick Hayden, Debbie W. Leung and Andreas Winter, Aspects of generic entanglement. 

Comm. Math. Phys. 265, 95-117 (2006). 
[21] Holger Dette, Strong approximation of eigenvalues of large dimensional Wishart matrices by roots 

of generalized Laguerre polynomials, J. Approx. Theory 118, 290-304 (2002). 



